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1. INTRODUCTION 

In this paper we report briefly the results of an extensive set of 
3-D hydrodynamic calculations that have been performed during the past 
two and one-half years to investigate the susceptibility of rotating 
clouds to gravitational fragmentation. Because of the immensity of 
parameter space and the expense of computations, we have chosen to 
restrict this investigation to strictly isothermal collapse sequences. 
Isothermality is, fortunately, a fairly accurate description of the 
optically thin phases of collapse of protostellar clouds: a one solar 
mass, Jeans unstable interstellar cloud, having a temperature of 10 K 
and density * 10® particles/cm 2 3 , for example, will collapse isothermal- 
ly through about 6 orders of magnitude in density before heating up 
significantly (Gerola and Glassgold, 1978). The assumption also ap- 
plies to protogalactic clouds (Silk, 1977). The results, therefore, 
should add specifically to our understanding of the dynamics of the 
earliest phases in a protostellar (or protogalactic) cloud's evolution. 
A much more complete discussion of these calculations can be found in 
Bodenheimer, Tohline and Black (1980), hereafter referred to as BTB. 


2. THE NUMERICAL CODE 

The 3-D hydrodynamic computer code that was used in this study is 
described in detail by Tohline (1980). Briefly, the code gives an 
explicit, first-order-accurate time-integration of the three-dimensional 
equation of motion and continuity equation on a moving-Eulerian, cylin- 
drical computational grid. A solution of the three-dimensional Poisson 
• equation at each time step permits the analysis of self-gravitating gas 
flows. In the model evolutions discussed here, a grid resolution of 
(32,16,32) in (R,6,Z) was used. The Z-axis was the rotation axis of 


each cloud. The 32 zones in the Z-direction resolved only the 
"northern" hemisphere of each cloud, as reflection symmetry through the 
equatorial plane was assusied in all cases. In most evolutions, the 16 
azimuthal zones spanned the entire 2n radians; in some cases, however, 
the 16 zones were distributed over only It radians and a periodic boun- 
dary condition was used in order to double the azimuthal resolution— of 
course, only "even" non-axisyumetric modes could be followed in these 
"n- symmetry" runs. No magnetic fielus or viscous forces are included 
in these calculations. 


3. THE CALCULATIONS 
3.1. Initial Models 

Each initial - model was ■- chosen to be a uniform-density (p ), 
-^uniform- temperature, spherical cloud in uniform rotation, on which a 
"well-defined non-axi symmetric density perturbation^ tfa> imposed. The 
■odels reported here were 1 M~ clouds of pure mbleCulpr fiy^^eo atTa 
temperature of 10 K; under the assumed isothermal behaviOrV itowevar, 
the results can be scaled to any mass and temperature desired. 

The axisymmetric initial models can be uniquely defined by a 
choice of the two parameters: 



total thermal energy 

(total gravitational potential energy | 


« _ total rotational kinetic energy 
** ~ | total gravitational potential energy | 


All of the discussion which follows centers .-around an analysis of 
whether or not gravitationally driven fragmentation in protostellar 
clouds is sensitive to variations in either one (or both) of these 
parameters . 

In all of the model evolutions reported here, the density pertur- 
bation imposed on the initial models had the following form: 


p * p Q (1.0 ♦ A 0 *f(R,Z)*cos(26)) 


( 1 ) 


The function f(R,Z), whose exact description is given in BTB, simply 
produced two orbiting Gauss ian-shaped "blobs" of maximum amplitude A. 
off-axis in the equatorial plane of the cloud. A. was chosen to be 
either 0.10 (10%) or 0.50 (50%) in these models. 


3.2. Typical (Axisynmetric) Evolution* 


Each aodel, being initially Jeans unstable to collapse, was fol- 
lowed for 1. 5-2.0 initial free-fall times (t f _). The calculations were 
stopped in each case when density gradients z $ecame to severe in a part 
of the computational grid- -hence, resolution became sufficiently poor — 
that a physically realistic representation of the fluid flow was no 
longer possible. Although these models include initial deviations from 
axial symmetry, little or no growth of oon-axisymmetric features took 
place during the fir^t free-fall time. As a result, the general be- 
havior of each evolution is similar to that described by axisynmetric 
(2-D) models (see, e.g., Black and Bodenheimer, 1976). During the 
first free-fall time, each model evolved through a sequence of succes- 
sively more flattened, centrally condensed spheroids, as rotational 
forces slowed collapse perpendicular to the rotation axis. At roughly 
1 1„, the cloud was sufficiently flat that isothermal pressure gradi- 
ents 1 were able to stop the collapse in the Z-direction, forming a 
centrally condensed disk structure. The meridional cross-section of 
one cloud, shown in figure 1, illustrates this structure. The axis 
ratio of the disk correlated with the initial value of a— a smaller a 
leading to a flatter disk— while the absolute size of the disk corre- 
lated with p. If forceo to remain axisynmetric, an off-axis density 
enhancement appeared in the disk and grew, via gravitational forces, to 
a well-defined toroidal (ring) structure with a density minimum at the 
cloud center. This ring grew in mass toward an "equilibrium” configu- 
ration analogous to that described by Ostriker (1964) for isothermal 
toruses. The ring, upon accretion of more material from the cloud, 
usually evolved beyond a critical mass-per-unit-leugth and collapsed 
upon itself toward an infinitesimally thin hoop. 


Figure 1. Density contours 
and velocity vectors in the 
meridional (R,Z) plane 
after 1.21 initial free- 
fall times for a model with 
a - 0.5 and p = 0.1. Maxi- 
mum _ density, 5 x 10 14 
g cm” 3 ; contour interval, 
factor of 3.2. Most of the 
cloud mass is in the flat- 
tened disk and the ring 
structure is evident. 




3.3. Non-Axi symmetric Behavior 


' £ 


The evolution of non-axisymmetric features in these isothermal 
clouds can be described quite simply. For high a (>.3) clouds— those 
clouds in which thermal pressure was initially important— the --- initial 
non-axisymmetric perturbation (NAP) underwent pressure- damping during 
the first free-fall time and the cloud evolved to a nearly axisymmetric 
ring structure. In all cases, the NAP began to grow after the ring 
formed, and in most instances this led to fragmentation of the ring 
into an equal mass binary system. A binary originating from "ring 
fragmentation" of this nature is shown in figure 2. In a few cases 
(a = .5), the ring began to collapse axisymmetrically toward a thin 
hoop before the NAP had had time to amplify significantly. In these 
instances, no fragmentation was evident at the time calculations were 
stopped. It should be emphasized however, that fragmentation of the 
ring was observed in all other high or models — models with a > .5 as 
well as those with a < .5. 

In low a models, the NAP did not damp during the initial collapse 
phase and in most cases actually amplified somewhat. Upon formation 
of the disk, the NAP evolved directly into a well-defined, equal mass 
binary system without going through an intermediate ring formation 
stage. A binary system orginating from this type of "blob" fragmen- 
tation is illustrated in figure 3. 




Figure 2 Figure 3 

Density contours in the equatorial plane after "ring” fragmentation has 
occurred in one model (Fig. 2) and after "blob" fragmentation has 
occurred in ^another model (Fig. 3). Points of maximum density 
(* 10” 13 g cm" 8 ) are marked by "M"; contour interval, factor of 1.7. 


Of the 35 complete evolutions that we have followed, 12 started 
from initial models that had a Gaussian-shaped NAP of amplitude 50%. 
The final state of these 12 models is summarised pictorially in figure 
4. A symbol is plotted for each model in this cr-0 plane, pinpointing 
its initial values of o and 0. Squares locate models whose final 
configuration was a well-defined binary system (the number below each 
symbol tells the density contrast between a fragment's maximum and its 
surrounding); circles locate models that did not fragment signifi- 
cantly. "Marginal stability" occurs around a = 0.35. This figure 
points out quite strongly two things: 


The degree to which fragmentation occurs, that is, the rate 
at which a NAP grows, is sensitive to a. 


Fragmentation of these isothermal clouds is practically 
independent of 0. 


Figure 4. The (a,0) plane 
indicating, for each initial 
a, 0, the final outcome of 
the evolution. Symbols are 
explained in the text. 


■Oft QAIMUAN KRTUMATION 



Several models were run from identical initial conditions but with 
"n- symmetry" imposed, in order to double the resolution in the azimu- 
thal coordinate. These evolutions showed that coarse zoning can lead 
to some artificial numerical damping of non-axisymmetric features, but 
that the qualitative results described above are correct. With im- 
proved resolution and, hence, less numerical damping, only a couple of 
minor modifications arose: (1) The "marginal stability" line in the 
a-0 plane moved to slightly higher o, and (2) models near the marginal 
stability line that had previously fragmented through the ring mode, 
changed and fragmented directly via the "blob" mode. 


Several models were also run with f(R,Z) = 1 (see eq. 1) and with 
a 10% initial NAP amplitude rather than the general 50% amplitude; all 
of these runs also used "n- symmetry" in order to avoid excessive 
damping. The lower amplitude models also evolved qualitatively in the 
manner described above. The only differences were minimal and were 
related to the fact that at a given time in an evolution, the amplitude 



of the NAP wee (understandably) lower. Fragmentation into a binary 
system was the general outcome. 


4. SUMMARY 

From the extensive set of numerical calculations briefly described 
above, it seems apparent that rotating, isothermal gas clouds are un- 
stable to fragmentation under a wide range of conditions. (CAUTION: 
This result for isothermal clouds cannot be generalized to all clouds, 
as is shown, for example, by Boss's analysis [these proceedings] of the 
stability of collapsing, adiabatic clouds.) It is of importance to 
note, however, that no fragmentation is apparent during a cloud's 
initial dynamic collapse toward a disk structure; rather it is the 
rotationally flattened disk/ring configuration that undergoes fragmen- 
tation. This is a considerably different picture ox fragmentation than 
has been presented, for example, by Hoyle (1953). 

The degree of instability and the mode (ring vs. blob) of fragmen- 
tation is sensitive to a, but insensitive to 0. The initial aoiplitude 
of a perturbation does not appear to be crucial— fragmentation should 
occur eventually even for low amplitude initial NAPs. 

Finally, it is of some interest to know what the properties are of 
the fragments that break out of these isothermal clouds. Before out- 
lining these properties we emphasize that in this set of calculations 
we have specifically excited the m * 2 (binary) non-axisymmetric mode; 
hence we have in some sense suppressed the development of other modes 
and we have promoted the development of equal mass components in the 
binary systems. In these evolutions, a typical fragment contained 
~ 15% of the initial cloud mass; had a specific angular momentum 
2 25-30% that of the original cloud; had a ratio of spin angular momen- 
tum to orbital angular momentum S 0.2; and itself had a ratio of ther- 
mal to gravitational energy or. <0.1. The formation of a binary 
system has therefore resulted in a 8 conversion of some of the original 
cloud's spin angular momentum into orbital angular momentum, and has 
produced "protostars" with reduced specific angular momenta . It is 
also evident that each fragment is unstable to further collapse (having 
low a) under the isothermal assumptions imposed here. 
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